Automatic Preprocessing Pipeline for White Matter Functional Analyses of Large-Scale Databases.

Recently, increasing evidence suggests that fMRI signals in white matter (WM), conventionally ignored as nuisance, are robustly detectable using appropriate processing methods and are related to neural activity, while changes in WM with aging and degeneration are also well documented. These findings suggest variations in patterns of BOLD signals in WM should be investigated. However, existing fMRI analysis tools, which were designed for processing gray matter signals, are not well suited for large-scale processing of WM signals in fMRI data. We developed an automatic pipeline for high-performance preprocessing of fMRI images with emphasis on quantifying changes in BOLD signals in WM in an aging population. At the image processing level, the pipeline integrated existing software modules with fine parameter tunings and modifications to better extract weaker WM signals. The preprocessing results primarily included whole-brain time-courses, functional connectivity, maps and tissue masks in a common space. At the job execution level, this pipeline exploited a local XNAT to store datasets and results, while using DAX tool to automatic distribute batch jobs that run on high-performance computing clusters. Through the pipeline, 5,034 fMRI/T1 scans were preprocessed. The intraclass correlation coefficient (ICC) of test-retest experiment based on the preprocessed data is 0.52 - 0.86 (N=1000), indicating a high reliability of our pipeline, comparable to previously reported ICC in gray matter experiments. This preprocessing pipeline highly facilitates our future analyses on WM functional alterations in aging and may be of benefit to a larger community interested in WM fMRI studies.

Functional alterations in bipartite network of white and grey matters during aging.

The effects of normal aging on functional connectivity (FC) within various brain networks of gray matter (GM) have been well-documented. However, the age effects on the networks of FC between white matter (WM) and GM, namely WM-GM FC, remains unclear. Evaluating crucial properties, such as global efficiency (GE), for a WM-GM FC network poses a challenge due to the absence of closed triangle paths which are essential for assessing network properties in traditional graph models. In this study, we propose a bipartite graph model to characterize the WM-GM FC network and quantify these challenging network properties. Leveraging this model, we assessed the WM-GM FC network properties at multiple scales across 1,462 cognitively normal subjects aged 22-96 years from three repositories (ADNI, BLSA and OASIS-3) and investigated the age effects on these properties throughout adulthood and during late adulthood (age ≥70 years). Our findings reveal that (1) heterogeneous alterations occurred in region-specific WM-GM FC over the adulthood and decline predominated during late adulthood; (2) the FC density of WM bundles engaged in memory, executive function and processing speed declined with age over adulthood, particularly in later years; and (3) the GE of attention, default, somatomotor, frontoparietal and limbic networks reduced with age over adulthood, and GE of visual network declined during late adulthood. These findings provide unpresented insights into multi-scale alterations in networks of WM-GM functional synchronizations during normal aging. Furthermore, our bipartite graph model offers an extendable framework for quantifying WM-engaged networks, which may contribute to a wide range of neuroscience research.

Comparing voxel- and feature-wise harmonization of complex graph measures from multiple sites for structural brain network investigation of aging.

Complex graph theory measures of brain structural connectomes derived from diffusion weighted images (DWI) provide insight into the network structure of the brain. Further, as the number of available DWI datasets grows, so does the ability to investigate associations in these measures with major biological factors, like age. However, one key hurdle that remains is the presence of scanner effects that can arise from different DWI datasets and confound multisite analyses. Two common approaches to correct these effects are voxel-wise and feature-wise harmonization. However, it is still unclear how to best leverage them for graph-theory analysis of an aging population. Thus, there is a need to better characterize the impact of each harmonization method and their ability to preserve age related features. We investigate this by characterizing four complex graph theory measures (modularity, characteristic path length, global efficiency, and betweenness centrality) in 48 participants aged 55 to 86 from Baltimore Longitudinal Study of Aging (BLSA) and Vanderbilt Memory and Aging Project (VMAP) before and after voxel- and feature-wise harmonization with the Null Space Deep Network (NSDN) and ComBat, respectively. First, we characterize across dataset coefficients of variation (CoV) and find the combination of NSDN and ComBat causes the greatest reduction in CoV followed by ComBat alone then NSDN alone. Second, we reproduce published associations of modularity with age after correcting for other covariates with linear models. We find that harmonization with ComBat or ComBat and NSDN together improves the significance of existing age effects, reduces model residuals, and qualitatively reduces separation between datasets. These results reinforce the efficiency of statistical harmonization on the feature-level with ComBat and suggest that harmonization on the voxel-level is synergistic but may have reduced effect after running through the multiple layers of the connectomics pipeline. Thus, we conclude that feature-wise harmonization improves statistical results, but the addition of biologically informed voxel-based harmonization offers further improvement.

Associations of baseline and longitudinal change in cerebellum volume with age-related changes in verbal learning and memory.

The cerebellum is involved in higher-order cognitive functions, e.g., learning and memory, and is susceptible to age-related atrophy. Yet, the cerebellum's role in age-related cognitive decline remains largely unknown. We investigated cross-sectional and longitudinal associations between cerebellar volume and verbal learning and memory. Linear mixed effects models and partial correlations were used to examine the relationship between changes in cerebellum volumes (total cerebellum, cerebellum white matter [WM], cerebellum hemisphere gray matter [GM], and cerebellum vermis subregions) and changes in verbal learning and memory performance among 549 Baltimore Longitudinal Study of Aging participants (2,292 visits). All models were adjusted by baseline demographic characteristics (age, sex, race, education), and APOE e4 carrier status. In examining associations between change with change, we tested an additional model that included either hippocampal (HC), cuneus, or postcentral gyrus (PoCG) volumes to assess whether cerebellar volumes were uniquely associated with verbal learning and memory. Cross-sectionally, the association of baseline cerebellum GM and WM with baseline verbal learning and memory was age-dependent, with the oldest individuals showing the strongest association between volume and performance. Baseline volume was not significantly associated with change in learning and memory. However, analysis of associations between change in volumes and changes in verbal learning and memory showed that greater declines in verbal memory were associated with greater volume loss in cerebellum white matter, and preserved GM volume in cerebellum vermis lobules VI-VII. The association between decline in verbal memory and decline in cerebellar WM volume remained after adjustment for HC, cuneus, and PoCG volume. Our findings highlight that associations between cerebellum volume and verbal learning and memory are age-dependent and regionally specific.

Joint independent component analysis for hypothesizing spatiotemporal relationships between longitudinal gray and white matter changes in preclinical Alzheimer's disease.

Characterizing relationships between gray matter (GM) and white matter (WM) in early Alzheimer's disease (AD) would improve understanding of how and when AD impacts the brain. However, modeling these relationships across brain regions and longitudinally remains a challenge. Thus, we propose extending joint independent component analysis (jICA) into spatiotemporal modeling of regional cortical thickness and WM bundle volumes leveraging multimodal MRI. We jointly characterize these GM and WM features in a normal aging (n=316) and an age- and sex-matched preclinical AD cohort (n=81) at each of two imaging sessions spaced three years apart, training on the normal aging population in cross-validation and interrogating the preclinical AD cohort. We find this joint model identifies reproducible, longitudinal changes in GM and WM between the two imaging sessions and that these changes are associated with preclinical AD and are plausible considering the literature. We compare this joint model to two focused models: (1) GM features at the first session and WM at the second and (2) vice versa. The joint model identifies components that correlate poorly with those from the focused models, suggesting the different models resolve different patterns. We find the strength of association with preclinical AD is improved in the GM to WM model, which supports the hypothesis that medial temporal and frontal thinning precedes volume loss in the uncinate fasciculus and inferior anterior-posterior association fibers. These results suggest that jICA effectively generates spatiotemporal hypotheses about GM and WM in preclinical AD, especially when specific intermodality relationships are considered a priori.

Herpes Viruses in the Baltimore Longitudinal Study of Aging: Associations With Brain Volumes, Cognitive Performance, and Plasma Biomarkers.

BACKGROUND AND OBJECTIVES: Although an infectious etiology of Alzheimer's Disease (AD) has received renewed attention with a particular focus on herpes viruses, the longitudinal effects of symptomatic herpes viruses (sHHV) infection on brain structure and cognition remain poorly understood, as does the effect of sHHV on AD/neurodegeneration biomarkers. METHODS: We used a longitudinal, community-based cohort to characterize the association of sHHV diagnoses with changes in 3T MRI brain volume and cognitive performance. Additionally, we related sHHV to cross-sectional differences in plasma biomarkers of AD (Aß42/40), astrogliosis (glial fibrillary acidic protein [GFAP]) and neurodegeneration (neurofilament light [NfL]). Baltimore Longitudinal Study of Aging (BLSA) participants were recruited from the community and assessed with serial brain MRIs and cognitive exams over an average of 3.4 (SD=3.2) and 8.6 (SD=7.7) years, respectively. sHHV classification used ICD9 codes documented at comprehensive health and functional screening evaluations at each study visit. Linear mixed effects and multivariable linear regression models were used in analyses. RESULTS: A total of 1,009 participants were included in the primary MRI analysis, 98% of whom were cognitively normal at baseline MRI (mean age = 65.7 years; 54.8% female). Having a sHHV diagnosis (N=119) was associated with longitudinal reductions in white matter volume (annual additional rate of change -0.34 cm3/year; p = 0.035), particularly in the temporal lobe. However, there was no association between sHHV and change in total brain, total gray matter, or AD signature region volume. Among the 119 participants with sHHV, exposure to antiviral treatment attenuated declines in occipital white matter (p = 0.04). Although the sHHV group had higher cognitive scores at baseline, sHHV diagnosis was associated with accelerated longitudinal declines in attention (annual additional rate of change -0.01 Z-score/year; p = 0.008). Additionally, sHHV diagnosis was associated with elevated plasma GFAP, but not related to Aß42/40 and NfL levels. DISCUSSION: These findings suggest an association of sHHV infection with white matter volume loss, attentional decline, and astrogliosis. Although the findings link sHHV to several neurocognitive features, the results do not support an association between sHHV and AD-specific disease processes.

Contrastive semi-supervised harmonization of single-shell to multi-shell diffusion MRI.

Diffusion weighted MRI (DW-MRI) harmonization is necessary for multi-site or multi-acquisition studies. Current statistical methods address the need to harmonize from one site to another, but do not simultaneously consider the use of multiple datasets which are comprised of multiple sites, acquisitions protocols, and age demographics. This work explores deep learning methods which can generalize across these variations through semi-supervised and unsupervised learning while also learning to estimate multi-shell data from single-shell data using the Multi-shell Diffusion MRI Harmonization Challenge (MUSHAC) and Baltimore Longitudinal Study on Aging (BLSA) datasets. We compare disentanglement harmonization models, which seek to encode anatomy and acquisition in separate latent spaces, and a CycleGAN harmonization model, which uses generative adversarial networks (GAN) to perform style transfer between sites, to the baseline preprocessing and to SHORE interpolation. We find that the disentanglement models achieve superior performance in harmonizing all data while at the same transforming the input data to a single target space across several diffusion metrics (fractional anisotropy, mean diffusivity, mean kurtosis, primary eigenvector).

Longitudinal changes of connectomes and graph theory measures in aging.

Changes in brain structure and connectivity in aging can be probed through diffusion weighted MRI and summarized with structural connectome matrices. Complex network analysis based on graph theory has been applied to provide measures that are correlated with neurobiological variations and can help guide quantitative study of brain function. However, the understanding of how connectomes change longitudinally is limited. In this work, we evaluate modern pipelines to obtain the connectomics data from diffusion weighted MRI scans across different sessions from control subjects (55-99 years old) in the Baltimore Longitudinal Study of Aging and Cambridge Centre for Ageing and Neuroscience. From the connectivity matrices, we compute graph theory measures to understand their brain networks and apply a linear mixed-effects model to study their longitudinal changes with respect to age. With this approach, we computed 14 graph theory measures of 1469 healthy subjects (2476 scans) and found statistically significant correlations between all 14 measures and age. In this analysis: 1) the brain becomes more segregated but less integrated in aging; 2) the overall network cost increases for older subjects; 3) the small-world organizations remain stable; and 4) due to high intra-subject variance, there is not clear trend for longitudinal changes of graph theory measures of individual subjects. Therefore, while useful to investigate brain evolution in aging at the population level, improvements in the connectome reconstruction are needed to decrease single subject variability for individual inference.

Aging and white matter microstructure and macrostructure: a longitudinal multi-site diffusion MRI study of 1218 participants.

Quantifying the microstructural and macrostructural geometrical features of the human brain's connections is necessary for understanding normal aging and disease. Here, we examine brain white matter diffusion magnetic resonance imaging data from one cross-sectional and two longitudinal data sets totaling in 1218 subjects and 2459 sessions of people aged 50-97 years. Data was drawn from well-established cohorts, including the Baltimore Longitudinal Study of Aging data set, Cambridge Centre for Ageing Neuroscience data set, and the Vanderbilt Memory & Aging Project. Quantifying 4 microstructural features and, for the first time, 11 macrostructure-based features of volume, area, and length across 120 white matter pathways, we apply linear mixed effect modeling to investigate changes in pathway-specific features over time, and document large age associations within white matter. Conventional diffusion tensor microstructure indices are the most age-sensitive measures, with positive age associations for diffusivities and negative age associations with anisotropies, with similar patterns observed across all pathways. Similarly, pathway shape measures also change with age, with negative age associations for most length, surface area, and volume-based features. A particularly novel finding of this study is that while trends were homogeneous throughout the brain for microstructure features, macrostructural features demonstrated heterogeneity across pathways, whereby several projection, thalamic, and commissural tracts exhibited more decline with age compared to association and limbic tracts. The findings from this large-scale study provide a comprehensive overview of the age-related decline in white matter and demonstrate that macrostructural features may be more sensitive to heterogeneous white matter decline. Therefore, leveraging macrostructural features may be useful for studying aging and could facilitate comparisons in a variety of diseases or abnormal conditions.

Accelerated decline in white matter microstructure in subsequently impaired older adults and its relationship with cognitive decline.

Little is known about a longitudinal decline in white matter microstructure and its associations with cognition in preclinical dementia. Longitudinal diffusion tensor imaging and neuropsychological testing were performed in 50 older adults who subsequently developed mild cognitive impairment or dementia (subsequently impaired) and 200 cognitively normal controls. Rates of white matter microstructural decline were compared between groups using voxel-wise linear mixed-effects models. Associations between change in white matter microstructure and cognition were examined. Subsequently impaired individuals had a faster decline in fractional anisotropy in the right inferior fronto-occipital fasciculus and bilateral splenium of the corpus callosum. A decline in right inferior fronto-occipital fasciculus fractional anisotropy was related to a decline in verbal memory, visuospatial ability, processing speed and mini-mental state examination. A decline in bilateral splenium fractional anisotropy was related to a decline in verbal fluency, processing speed and mini-mental state examination. Accelerated regional white matter microstructural decline is evident during the preclinical phase of mild cognitive impairment/dementia and is related to domain-specific cognitive decline.

Recognition Memory is Associated with Distinct Patterns of Regional Gray Matter Volumes in Young and Aged Monkeys.

Cognitive aging varies tremendously across individuals and is often accompanied by regionally specific reductions in gray matter (GM) volume, even in the absence of disease. Rhesus monkeys provide a primate model unconfounded by advanced neurodegenerative disease, and the current study used a recognition memory test (delayed non-matching to sample; DNMS) in conjunction with structural imaging and voxel-based morphometry (VBM) to characterize age-related differences in GM volume and brain-behavior relationships. Consistent with expectations from a long history of neuropsychological research, DNMS performance in young animals prominently correlated with the volume of multiple structures in the medial temporal lobe memory system. Less anticipated correlations were also observed in the cingulate and cerebellum. In aged monkeys, significant volumetric correlations with DNMS performance were largely restricted to the prefrontal cortex and striatum. Importantly, interaction effects in an omnibus analysis directly confirmed that the associations between volume and task performance in the MTL and prefrontal cortex are age-dependent. These results demonstrate that the regional distribution of GM volumes coupled with DNMS performance changes across the lifespan, consistent with the perspective that the aged primate brain retains a substantial capacity for structural reorganization.

Disease Burden Affects Aging Brain Function.

BACKGROUND: Most older adults live with multiple chronic disease conditions, yet the effect of multiple diseases on brain function remains unclear. METHODS: We examine the relationship between disease multimorbidity and brain activity using regional cerebral blood flow (rCBF) 15O-water PET scans from 97 cognitively normal participants (mean baseline age 76.5) in the Baltimore Longitudinal Study of Aging (BLSA). Multimorbidity index scores, generated from the presence of 13 health conditions, were correlated with PET data at baseline and in longitudinal change (n = 74) over 5.05 (2.74 SD) years. RESULTS: At baseline, voxel-based analysis showed that higher multimorbidity scores were associated with lower relative activity in orbitofrontal, superior frontal, temporal pole and parahippocampal regions, and greater activity in lateral temporal, occipital, and cerebellar regions. Examination of the individual health conditions comprising the index score showed hypertension and chronic kidney disease individually contributed to the overall multimorbidity pattern of altered activity. Longitudinally, both increases and decreases in activity were seen in relation to increasing multimorbidity over time. These associations were identified in orbitofrontal, lateral temporal, brainstem, and cerebellar areas. CONCLUSION: Together, these results show that greater multimorbidity is associated with widespread areas of altered brain activity, supporting a link between health and changes in aging brain function.

Validation of Group-wise Registration for Surface-based Functional MRI Analysis.

Resting-state functional MRI (rsfMRI) provides important information for studying and mapping the activities and functions of the brain. Conventionally, rsfMRIs are often registered to structural images in the Euclidean space without considering cortical surface geometry. Meanwhile, a surface-based representation offers a relaxed coordinate chart, but this still requires surface registration for group-wise data analysis. In this work, we investigate the performance of two existing surface registration methods in a surface-based rsfMRI analysis framework: FreeSurfer and Hierarchical Spherical Deformation (HSD). To minimize registration bias, we establish shape correspondence using both methods in a group-wise manner that estimates the unbiased average of a given cohort. To evaluate their performance, we focus on neuroanatomical alignment as well as the amount of distortion that can potentially bias surface tessellation for secondary level rsfMRI data analyses. In the pilot analysis, we examine a single timepoint of imaging data from 100 subjects out of an aging cohort. Overall, HSD establishes improved shape correspondence with reduced mean curvature deviation (10.94% less on average per subject, paired t-test: p <10-10) and reduced registration distortion (FreeSurfer: average 41.91% distortion per subject, HSD: 18.63%, paired t-test: p <10-10). Furthermore, HSD introduces less distortion than FreeSurfer in the areas identified in the individual components that were extracted by surface-based independent component analysis (ICA) after spatial smoothing and time series normalization. Consequently, we show that FreeSurfer capture individual components with globally similar but locally different patterns in ICA in visual inspection.

Joint cortical surface and structural connectivity analysis of Alzheimer's Disease.

Prior neuroimaging studies have demonstrated isolated structural and connectivity changes in the brain due to Alzheimer's Disease (AD). However, how these changes relate to each other is not well understood. We present a preliminary study to begin to fill this gap by leveraging joint independent component analysis (jICA). We explore how jICA performs in an analysis of T1 and diffusion weighted MRI by characterizing the joint changes of complex cortical surface and structural connectivity metrics in AD in subjects from the Baltimore Longitudinal Study of Aging. We calculate 588 region-based cortical metrics and 4,753 fractional anisotropy-based connectivity metrics and project them into a low-dimensional manifold with principal component analysis. We perform jICA on the manifold and subsequently backproject the independent components to the original data space. We demonstrate component stability with 3-fold cross validation and find differential component loadings between 776 cognitively unimpaired control subjects and 23 with AD that generalizes across folds. In addition, we perform the same analysis on the surface and connectivity metrics separately and find that the joint approach identifies both novel and similar components to the separate approaches. To illustrate the joint approach's primary utility, we provide an example hypothesis for how surface and connectivity components may vary together with AD. These preliminary results suggest jointly varying independent cortical surface and structural connectivity components can be consistently extracted from MRI data and provide a data-driven way for generating novel hypotheses about AD that may not be captured by separate analyses.

Proof-of-concept evidence for trimodal simultaneous investigation of human brain function.

The link between spatial (where) and temporal (when) aspects of the neural correlates of most psychological phenomena is not clear. Elucidation of this relation, which is crucial to fully understand human brain function, requires integration across multiple brain imaging modalities and cognitive tasks that reliably modulate the engagement of the brain systems of interest. By overcoming the methodological challenges posed by simultaneous recordings, the present report provides proof-of-concept evidence for a novel approach using three complementary imaging modalities: functional magnetic resonance imaging (fMRI), event-related potentials (ERPs), and event-related optical signals (EROS). Using the emotional oddball task, a paradigm that taps into both cognitive and affective aspects of processing, we show the feasibility of capturing converging and complementary measures of brain function that are not currently attainable using traditional unimodal or other multimodal approaches. This opens up unprecedented possibilities to clarify spatiotemporal integration of brain function.

Default mode network connectivity and cognition in the aging brain: the effects of age, sex, and APOE genotype.

The default mode network (DMN) overlaps with regions showing early Alzheimer's Disease (AD) pathology. Age, sex, and apolipoprotein E ɛ4 are the predominant risk factors for developing AD. How these risk factors interact to influence DMN connectivity and connectivity-cognition relationships before the onset of impairment remains unknown. Here, we examined these issues in 475 cognitively normal adults, targeting total DMN connectivity, its anticorrelated network (acDMN), and the DMN-hippocampal component. There were four main findings. First, in the ɛ3 homozygous group, lower DMN and acDMN connectivity was observed with age. Second, sex and ɛ4 modified the relationship between age and connectivity for the DMN and hippocampus with ɛ4 vs. ɛ3 males showing sustained or higher connectivity with age. Third, in the ɛ3 group, age and sex modified connectivity-cognition relationships with the oldest participants having the most differential patterns due to sex. Fourth, ɛ4 carriers with lower connectivity had poorer cognitive performance. Taken together, our results show the three predominant risk factors for AD interact to influence brain function and function-cognition relationships.

PreQual: An automated pipeline for integrated preprocessing and quality assurance of diffusion weighted MRI images.

PURPOSE: Diffusion weighted MRI imaging (DWI) is often subject to low signal-to-noise ratios (SNRs) and artifacts. Recent work has produced software tools that can correct individual problems, but these tools have not been combined with each other and with quality assurance (QA). A single integrated pipeline is proposed to perform DWI preprocessing with a spectrum of tools and produce an intuitive QA document. METHODS: The proposed pipeline, built around the FSL, MRTrix3, and ANTs software packages, performs DWI denoising; inter-scan intensity normalization; susceptibility-, eddy current-, and motion-induced artifact correction; and slice-wise signal drop-out imputation. To perform QA on the raw and preprocessed data and each preprocessing operation, the pipeline documents qualitative visualizations, quantitative plots, gradient verifications, and tensor goodness-of-fit and fractional anisotropy analyses. RESULTS: Raw DWI data were preprocessed and quality checked with the proposed pipeline and demonstrated improved SNRs; physiologic intensity ratios; corrected susceptibility-, eddy current-, and motion-induced artifacts; imputed signal-lost slices; and improved tensor fits. The pipeline identified incorrect gradient configurations and file-type conversion errors and was shown to be effective on externally available datasets. CONCLUSIONS: The proposed pipeline is a single integrated pipeline that combines established diffusion preprocessing tools from major MRI-focused software packages with intuitive QA.

Hippocampal activation and connectivity in the aging brain.

The hippocampus and underlying cortices are highly susceptible to pathologic change with increasing age. Using an associative face-scene (Face-Place) encoding task designed to target these regions, we investigated activation and connectivity patterns in cognitively normal older adults. Functional MRI scans were collected in 210 older participants (mean age = 76.4 yrs) in the Baltimore Longitudinal Study of Aging (BLSA). Brain activation patterns were examined during encoding of novel Face-Place pairs. Functional connectivity of the hippocampus was also examined during encoding, with seed regions placed along the longitudinal axis in the head, body and tail of the structure. In the temporal lobe, task activation patterns included coverage of the hippocampus and underlying ventral temporal cortices. Extensive activation was also seen in frontal, parietal and occipital lobes of the brain. Functional connectivity analyses during overall encoding showed that the head of the hippocampus was connected to frontal and anterior/middle temporal regions, the body with frontal, widespread temporal and occipital regions, and the tail with posterior temporal and occipital cortical regions. Connectivity limited to encoding of subsequently remembered stimuli showed a similar pattern for the hippocampal body, but differing patterns for the head and tail regions. These results show that the Face-Place task produces activation along the occipitotemporal visual pathway including medial temporal areas. The connectivity results also show that patterns of functional connectivity vary throughout the anterior-posterior extent of the hippocampus during memory encoding. As these patterns include regions vulnerable to pathologic change in early stages of Alzheimer's disease, continued longitudinal assessment of these individuals can provide valuable information regarding changes in brain-behavior relationships that may occur with advancing age and the onset of cognitive decline.

Pandora: 4-D White Matter Bundle Population-Based Atlases Derived from Diffusion MRI Fiber Tractography.

Brain atlases have proven to be valuable neuroscience tools for localizing regions of interest and performing statistical inferences on populations. Although many human brain atlases exist, most do not contain information about white matter structures, often neglecting them completely or labelling all white matter as a single homogenous substrate. While few white matter atlases do exist based on diffusion MRI fiber tractography, they are often limited to descriptions of white matter as spatially separate "regions" rather than as white matter "bundles" or fascicles, which are well-known to overlap throughout the brain. Additional limitations include small sample sizes, few white matter pathways, and the use of outdated diffusion models and techniques. Here, we present a new population-based collection of white matter atlases represented in both volumetric and surface coordinates in a standard space. These atlases are based on 2443 subjects, and include 216 white matter bundles derived from 6 different automated state-of-the-art tractography techniques. This atlas is freely available and will be a useful resource for parcellation and segmentation.

Distortion correction of diffusion weighted MRI without reverse phase-encoding scans or field-maps.

Diffusion magnetic resonance images may suffer from geometric distortions due to susceptibility induced off resonance fields, which cause geometric mismatch with anatomical images and ultimately affect subsequent quantification of microstructural or connectivity indices. State-of-the art diffusion distortion correction methods typically require data acquired with reverse phase encoding directions, resulting in varying magnitudes and orientations of distortion, which allow estimation of an undistorted volume. Alternatively, additional field maps acquisitions can be used along with sequence information to determine warping fields. However, not all imaging protocols include these additional scans and cannot take advantage of state-of-the art distortion correction. To avoid additional acquisitions, structural MRI (undistorted scans) can be used as registration targets for intensity driven correction. In this study, we aim to (1) enable susceptibility distortion correction with historical and/or limited diffusion datasets that do not include specific sequences for distortion correction and (2) avoid the computationally intensive registration procedure typically required for distortion correction using structural scans. To achieve these aims, we use deep learning (3D U-nets) to synthesize an undistorted b0 image that matches geometry of structural T1w images and intensity contrasts from diffusion images. Importantly, the training dataset is heterogenous, consisting of varying acquisitions of both structural and diffusion. We apply our approach to a withheld test set and show that distortions are successfully corrected after processing. We quantitatively evaluate the proposed distortion correction and intensity-based registration against state-of-the-art distortion correction (FSL topup). The results illustrate that the proposed pipeline results in b0 images that are geometrically similar to non-distorted structural images, and more closely match state-of-the-art correction with additional acquisitions. In addition, we show generalizability of the proposed approach to datasets that were not in the original training / validation / testing datasets. These datasets included varying populations, contrasts, resolutions, and magnitudes and orientations of distortion and show efficacious distortion correction. The method is available as a Singularity container, source code, and an executable trained model to facilitate evaluation.